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Abstract 

We study the transport properties of a quantum dot driven by either a rotating magnetic field or 
an ac gate voltage using the Floquet master-equation approach. Both types of ac driving lead to 
photon-assisted tunneling where quantized amounts of energy are exchanged with the driving field. 
It is found that the differential-conductance peak due to photon-assisted tunneling does not survive 
in the Coulomb-blockade regime when the dot is driven by a rotating magnetic field. Furthermore, 
we employ a generalized MacDonald formula to calculate the time-averaged noise spectra of ac- 
driven quantum dots. Besides the peak at zero frequency, the noise spectra show additional peaks 
or dips in the presence of an ac field. For the case of an applied ac gate voltage, the peak or dip 
position is fixed at the driving frequency, whereas the position changes with increasing amplitude 
for the case of a rotating magnetic field. Additional features appear in the noise spectra if a dc 
magnetic field is applied in addition to a rotating field. In all cases, the peak or dip positions can 
be understood from the energy differences of two available Floquet channels. 

PACS numbers: 73.63.Kv, 73.23.Hk, 72.10.Bg, 05.60. Gg 
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I. INTRODUCTION 



Quantum conductors based on single molecules or semiconductor quantum dots are 
promising building blocks for future electronics and model systems for the study of funda- 
mental quantum phenomena. 1 However, much information on quantum conductors is beyond 
the reach of measurements of the current or conductance alone. Instead, the understanding 
of the transport properties calls for a study of the full counting statistics. 2-11 In the past 
decade, valuable information on microscopic details of the charge transport has been ob- 
tained from measurements of the current fluctuations or current noise. 12 Previous studies 
have shown that one can extract parameters such as the average backscattered charge, 13 the 
intrinsic time scales, 14,15 and the asymmetry of the dot-lead coupling 16,17 from current-noise 
measurements. Most studies of the current noise have focused on the zero-frequency noise 
power S(0). 18-20 The zero- frequency noise reflects the average properties of the tunneling. 
Since the finite-frequency current noise S(u) is a measure of the correlations between tun- 
neling events with their time difference conjugate to the frequency u 21 it is interesting to go 
beyond the zero-frequency limit. Aguado and Brandes 14 have demonstrated that the noise 
spectra can show dip structures at the splitting energy of an open quantum two-level sys- 
tem, with their width controlled by its dissipative dynamics. In many works, the MacDonald 
formula 22 has been used to study the current-noise spectra. 23-25 

Effective in-situ manipulation of quantum conductors is a key step for further develop- 
ment. Using a time- dependent field to manipulate the dynamics of quantum dots promises 
to be advantageous in situations ranging from photon-assisted inelastic tunneling 26 to quan- 
tum pumping. 27 When the conductor is driven by an ac field, one expects novel features 
due to the interplay of intrinsic oscillation frequencies and the external driving frequency. 
Several recent studies 28,29 indicate that key information is hidden in the noise spectra of 
the ac-driven transport. For instance, Barrett and Stace 28 have proposed to extract the 
characteristic timescales such as the inverse dephasing and relaxation rates of a solid-state 
charge qubit coupled to a microwave field from the noise spectrum. Wabnig et al. 29 have 
proposed to estimate the coherence time of the spin in a quantum dot by measuring its noise 
spectra under an ac magnetic field. These results are obtained based on the rotating-wave 
approximation and usually in the limit of infinite on-site Coulomb interaction. 

For periodically driven systems, an appropriate theoretical tool to go beyond the rotating- 
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wave approximation is the Floquet theorem. 30,31 Various attempts have been made by gen- 
eralizing the existing steady-state transport approaches such as the scattering matrix 32-34 
and nonequilibrium Green's functions 35-37 with the help of the Floquet theorem. However, 
these methods are not adequate to fully take the Coulomb blockade in quantum dots into 
account, which dominates the transport properties of small-size quantum conductors. The 
quantum master equation 38 in its various manifestations 39,40 is able to give a good account 
of the Coulomb blockade in the weak-tunneling limit. This method has previously been 
generalized using the Floquet theorem to study the current and the zero-frequency noise 
power in an ac-driven conductor. 38,41-45 

In the present study, we employ the Floquet master equation in the Fock space of an ac- 
driven quantum dot to study the transport properties such as the differential conductance 
and the full noise spectrum in the sequential-tunneling limit. As the ac field we consider a 
rotating; magnetic field as well HS fill cLC gate voltage for comparison. We employ a generalized 
MacDonald formula for the time-averaged noise spectra in the presence of a periodic ac 
field. An equivalent form of the generalized MacDonald formula has been given by Clerk 
and Girvin 46 without derivation. For completeness, we present a derivation in the appendix. 
Note that the authors are concerned with a different case, namely an ac bias voltage, and 
do not employ the Floquet formalism. 

Our paper is organized as follows. In Sec. II, a Floquet master-equation formalism is 
presented to study the transport properties of ac-driven quantum dots. Expressions for the 
noise spectra are derived based on the full counting statistics and the generalized MacDonald 
formula. In Sec. Ill, the transport properties of the ac-driven quantum dot are studied. The 
ac field is either a rotating magnetic field or an oscillating on-site energy due to a periodic 
gate voltage. The characteristic features in the transport properties are presented and 
discussed. In Sec. IV, a brief summary is given. 

II. FORMALISM 
A. Model 

In this paper, we study the transport properties of a single-level quantum dot driven by 
an ac field. The quantum dot is coupled to the left (L) and right (R) electron leads. The 
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leads are assumed to be ideal and free of interactions. The Hamiltonian of the model system 
can be written as 

H(t) =H L + H R + H dot (t) + H T = H + H T , (1) 

where H& ot is the Hamiltonian of the isolated quantum dot, which contains the effects of 
the ac field and the Coulomb interaction, Hi = X)ka e «kcrcJ kcr Qkcr represents the Hamiltonian 
of lead I = L,R, where Qko- ( c Lo-) annihilates (creates) an electron with spin a, crystal 
momentum k, and energy e/ko- in lead I, and Ht = J2ika Vivurcj^do- + h.c. describes the 
coupling between the quantum dot and the leads, where cft a (d a ) is the spin-cr electron 
creation (annihilation) operator in the quantum dot. We note that we do not assume an 
infinite Coulomb-interaction strength U — > oo, in contrast to previous studies. 44 Instead, 
finite Coulomb interaction will be included by taking the doubly occupied state into account. 

In the following, we focus on the limit of weak dot-lead coupling and investigate the 
transport properties of the quantum dot using the Floquet master-equation method. For a 
small quantum dot, for which the Coulomb interaction can dominate the transport behavior, 
the treatment of the Coulomb interaction must go beyond the mean-field level. To this end, 
it is convenient to rewrite the dot Hamiltonian in the electron-number basis of the Fock 
space. 47 ' 48 In this description, the quantum dot can either be in the empty state |0), the 
singly occupied state \a) with spin a =f or J,, or the doubly occupied state | "fi.). In the 
following, we denote the states in the Fock space by Latin letters, \a) = |0), | f), | tl}- 
Using this orthonormal basis, the dot-lead coupling can be described naturally with the help 
of Hubbard operators 

X ab = \a)(b\, (2) 

which describes the transition of the quantum dot from state \b) to state \a). The second- 
quantized dot-electron creation operator can thus be rewritten in terms of the Hubbard 
operators as 

4 = \a){0\ + rj ff \UWl (3) 

where a represents the opposite spin of a and the factor r] a = ±1 for a =t, 4, respectively, is 
due to the anticommutation relation of the Fermions. In terms of these Hubbard operators, 
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the Hamiltonian for the dot-lead coupling and the isolated dot can be rewritten as 

H Aot = Y.H°\a)(bl (4) 

ab 

Ht= E V^4A + h.c., (5) 

ab,ik<7 

respectively. Here, we have made the coupling strength V$? c depend on the occupancy of 
the initial and final states of the transition in the quantum dot. The explicit form of the 
dot Hamiltonian depends on the details of the device geometry and the external ac driving 
field. It will be specified in the following sections. 



B. The Floquet quantum-master-equation approach 

1. Floquet states 

Due to the presence of a time-periodic external field, the dynamics of the quantum dot 
is governed by a Hamiltonian that is periodic in time with the frequency Q = 2tt/T, i.e., 
Hdot(t) = Hdot(t + T), where T denotes the period. The solution of the time-periodic 
Hamiltonian can be simplified by the Floquet theorem, 30,31 which states that the solution 
of the Schrodinger equation for the dot Hamiltonian can be obtained from (we set % — |e| = 
ks = 1 in the following) 

a(t)) =e a \a(t)), (6) 

where e a is the time-independent Floquet quasienergy and \a(t)) is the corresponding Flo- 
quet state, which has the same period T, \a(t)) = \a(t + T)). Here, Greek letters are used 
to denote the Floquet states. Further simplification is possible by decomposing the Floquet 
states into a Fourier series, 

|a(t)>=£e-^K), (7) 
k 

with the reverse transformation 

\a k ) = l- f dte lknt \a(t)) (8) 

and analogously for H& ot (t). The Fourier transform of Eq. (6) then reads 

£#dot,fc-fc' \otk') ~ kQ \a k ) = e a \a k ). (9) 
k' 
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H do t(t) 



The quasienergy e a can evidently be restricted to the first Brillouin zone [0, Q) of the Floquet 
space, while the Floquet index k can assume any integer value. Equivalently, we can view 
e a + kVt as the quasienergy in the extended zone scheme. 

We also introduce the Hubbard operator in the Floquet states to describe the transition 
between the Floquet states as X a p(t) = \a(t)){(3(t)\. For the time-dependent transport, 
it is more convenient to work with these Floquet states. This is most advantageous in 
transformations of the following form, which we will use in the derivation below, 

X a p(t',t) = U^t',t)X a p(t')U (t',t) 

= e^e-toW-Vx^t), (10) 

where 

U (t',t) = T c exp i^-ij' dt" [H L +H R + H dot (t")]\ (11) 

denotes the time-evolution operator due to the Hamiltonian in the absence of tunneling. 
Here, T c is the time-ordering operator and the dot Hamiltonian is explicitly time-dependent. 



2. The Floquet quantum master equation with counting fields 

For a quantum dot coupled to external leads, the exact quantum master equation can be 
written in the interaction picture as 38 ' 39 

j tPl {t) = -i[H T>I (t), Pl (t )] - jyt'[H T>I (t),[H T!l (t'), Pl (t')}}, (12) 

where Aj(t) = UQ(t,t )A(t)Uo(t,t Q ) denotes an operator in the interaction picture and pit) 
is the density matrix in the Fock space of the full system. 

A complete description of the electronic transport through the quantum dot is pro- 
vided by the full counting statistics. Properties such as the noise spectrum are deter- 
mined by the counting statistics of the electrons arriving at and departing from the leads. 
All information on the counting statistics is contained in the moment-generating function 
4>(Xl,Xr) — ( ex P(^XL^L + ^XrXr)) ■ Here, \i represents the counting field in the lead Z, 
which counts how many electrons have tunneled into or out of the lead. Ni = c aur c nur 
is the electron-number operator in lead /. We introduce the operator 

F(Xl, XR, t) = Tr lcads e^ N ^ N «p(t). (13) 
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In this we follow Kaiser and Kohler, 44 except that we introduce two counting fields. In the 
limit of xl — and \R ~ * 0, T becomes the reduced density matrix of the quantum dot, 
p dot = Tr i cads p. Moreover, the moment-generating function <f>(xL,XRjt) can be obtained by 
tracing out the dot degrees of freedom, <fi = Tr^ot J 7 - We decompose T into a Taylor series, 



' (lXL) m (tXRT 



T=yy y ALI y ^ KI j" m ' n , (14) 

m =0n=0 HI. lb. 

where the coefficients 

f)m+n 

= Tr leads N™N n R p (15) 



3m+n 

-pm,n _ jp 

' d m (t X L)d n (i XR ) 



provide a direct access to the moments (N^N^) = Tr^ot J zm ' n . In particular, we obtain the 
reduced density matrix of the quantum dot, pdot = J 70 ' - 

To find the solutions for J 7 , we first transform the equation of motion for the density 
matrix in the interaction picture, Eq. (12), back to the Schrodinger picture, 

^M + i[H (t),p(t)] = -i[H T ,tf (t ,t)p(t )U (t Q ,t)} 

- f 'dt'[H T ,tf (t',t)[H T ,p(t')}U (t',t)}. (16) 

Jt 

Then, we multiply by e lXLNL+lXRNR from the left and take the trace over the lead degrees of 
freedom to obtain 

— + i[H doi (t), jF{xl, Xr, t)\ 

= -I Tr lcads ^ L N L+ ixnN R ^ jjl^ t) p (t )U (to, t)] 

- f dt' Tr lcads jXLN L -HxHtf*\H^ C /t^ ) t)[H ^ p (t')]U (t', t)\, (17) 

which is still exact. 

We now assume that the full density operator is of product form at the initial time to, 
p(to) = pdot(^o) ® P? C ads5 where Pi cads describes the leads in separate thermal equilibrium. 
This assumption is reasonable since we are not interested in transient effects coming from 
the initial state. Such effects have been studied by Flindt et a/. 49 The first term on the 
right-hand side of Eq. (17) then vanishes. Furthermore, we make the sequential-tunneling 
approximation appropriate for weak tunneling, i.e., we treat the tunneling perturbatively to 
second order in H?- Since two powers of are already explicit in the second term on the 
right-hand side of Eq. (17), we can express p(t') in terms of the unperturbed time evolution, 
p{t') ~ Uo(t,t')p(t)U (t,t'). This makes the master equation local in time, i.e., Markovian. 



For details, see, e.g., Ref. 39. We thus do not include non-Markovian effects as studied by 
Flindt et a/. 49 This is valid if the relaxation time in the leads is short compared to the typical 
timescales of the dot, which in our case include the period of the ac field. Since relaxation 
times in metallic leads are of the order of femtoseconds, this is easily satisfied. 
Finally, we send t — oo and obtain 

dHXL dt XR,t) + ^dot(*)^(Xz, Xr, t)} (18) 

roo „ 

= - / drTr lead8 e M+i ^ NR [H T ,[H T (t-r,t),p(t)}}, 
Jo 

where H T {t' ,t) = U j (t',t)H T U (t',t). 

To obtain the Floquet master equation, we write Eq. (18) in the basis of Floquet states 
\a(t)), \j3(t)). By making use of the relation Eq. (10) and tracing out the lead degrees of 
freedom, we arrive at the equation of motion for J 7 , 

J^iXL, XR, t) = { [C + (e** - 1)J L+ + (e-** - 1)J L _ 

+ (e** - 1)J R+ + (e-'x* - 1)J R J\ F(XL,XR,t)} , (19) 

where the superoperators are given by 

(2 

in 7<5 l=L,R;cr 

(fi(e)e ier T[l(e) + Re^r 1 *) e^~^ ■<a(f)|a)<&| 7 (f)><7(f - r)\m)(n\8(t - r))^( X L,XR,t) 



(CJ 7 )^ = -i{e a - ep)T a p + — j dr J del 

^ ^ I ab;mn 7<5 l=L,R;cr 



+(/^)e i£r r^. mn ( £ ) + f^e'^T^.J e^-^>(t)|a)(6| 7 (t)> ^ 5 (xl, Xr, t)(6(t - r)\m)(n\f3(t - r)> 
+(Me)e- i£T T% mn (e) + h{e)e^ n . ah ) <X*—*{ a {t - r)|o)(6| 7 (t - r))^ s (x L , X R,t)(5(t)\m)(n\P(t)) 
-{Ue)e^Tt n , ab {e) + fi&e^H^) ?«r(XL, XR,t)(l(t - r)\a)(b\8(t - r))(5(t)\m)(n\f3(t)) 

roo p _ 

{Jl + FU = ^J drj deJ2 U^Xmn 
^ ^ ab;mn -yS cr 

(e ie V^- e ^(a(t)|a)(6| 7 (t))J- 75 ( X L,X R ,t)(5(t-r)|m)(n|/3(t-r)) 
+e-^e^-^(a(t - r)|a)(6| 7 (* - r))^ s ( X L,XR,t)(5(t)\m)(n\/3(t))) , (21) 

roo r 

ah Jo j u x 

ao;mn 70 & 

( e --V^-^>(t)|a>(^ 

+e l£ V^- £ ^(a(t-r)|a)(6| 7 (t-r))J- 7 ,( Xi , X/? ,t)(5(t)|m)(n|/3(t))). (22) 



Here, we have defined the tunneling rate T^ n . afe (e) = 2iipi(e)V l a ^ a V{l n * , where pi(e) is the 
density of states in lead I. 

Inserting the Taylor expansion of J 7 in Eq. (14) into its equation of motion [Eq. (19)], 



one obtains a hierarchy of equations for the expansion coefficients, 

= CF°>\ (23) 

jF 1 ' = CI** + (J L+ - J L ^ > , (24) 

J?'* 1 = C^ 1 + (J R+ - Jr-)^ , (25) 

= + 2{J L+ - Jl-)^ 1 ' + (J L+ + J L S)^\ (26) 

jT^ = £F°> 2 + 2(J R+ - J R _)^ 1 + (J R+ + J R -)F°'°, (27) 

j t T x > x = LT X > X + (J L+ - Jl-)^ ' 1 + (J R+ - J R -)F lfi , etc. (28) 

As described above, these coefficients contain the full counting statistics. The charge current 
out of lead / is defined as the negative of the time-derivative of the charge in lead I, Ii(t) = 
edNi/dt. The final expression for the current out of the left lead is given by 44 

(I L (t)} = e Tr dot (P>°) = e Tr dot ( J L+ - J L J) F°>°. (29) 



The dc component of the current then gives the time average I. In order to find J 70 ' , 
we have to solve a set of linear equations with the help of the normalization condition of 
probability Tr^ot J 70 ' ^) = 1. 

3. Generalized MacDonald formula for time-averaged noise spectra 

We are interested in the frequency- dependent current noise of the quantum dot driven 
by an ac field. The zero-frequency current noise for non-adiabatical driving has been in- 
vestigated in Ref. 44 using the Floquet master-equation approach in the Coulomb-blockade 
regime. The symmetrized current-current correlation function is defined by 

s lv (t,t') = (iimit')) + (iv(t')h(t)) - 2(i l (t))(i l ,(t')), (30) 

where Ii(t) represents the current operator at the time t from the lead /. The current-noise 
spectra are defined as the Fourier transform of S^/(t, t'). Since our system is driven by an ac 
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field, the current noise is a double-time function. However, the periodicity of our problem 
makes it possible to characterize the spectra by averaging over one driving period. 

At finite frequencies, the total current I(t) measured by a measurement device depends 
on both the particle and the displacement currents in the lead-dot-lead junction. If one 
expresses the displacement currents by the particle currents II and I R , one obtains the 
Ramo-Shockley theorem, 12 ' 50 ' 51 

I{t)=aI L (t)-bI R {t). (31) 

Here the coefficients a and b, which satisfy a + b = 1, are specified by the device geometry. 
It is straightforward to show that the total time-averaged noise spectrum is given by 

S(u>) = a 2 S LL (u) + b 2 S RR (co) - ab(S LR (u) + S RL (u)), (32) 

where Sw(oj) (1,1' = L,R) represents the frequency-dependent time-averaged current corre- 
lation between I t and ///, 

S w (u) = \\l J dt'e^-^SuitX). (33) 

In this study, we used two counting fields to derive the noise spectra. An alternative ap- 
proach is to calculate the charge fluctuation on the dot employing the quantum regression 
formula. 3,14 ' 52,53 The two approaches are physically equivalent due to the charge conservation 
condition in the transport. 

The formula for the zero-frequency noise has been presented in Ref. 44. For the two- 
terminal device, it is adequate to find the time-averaged zero-frequency noise from the 
fluctuations of the current flowing out of a chosen lead, 5(0) = Su(0). The solution for 5(0) 
resulting from the Floquet quantum master equation reads 

5(0) = % f T dt e 2 Tr dot \2(J l+ - Ji-)^ lAr + (J l+ + Ji-)F°'°] , (34) 
I Jo J 

where the prefactor of 2 is inserted to make the noise formula consistent with Ref. 12. For 

Poissonian noise, we then obtain 5(0) = 2el. 5ij is the usual Kronecker symbol. Following 

Ref. 44, the new function j^^- L < SlR in the noise expression is defined as 

jrWi* = jAlAr _ jro.o Tr dot {F SlL ' SlR }, (35) 

and satisfies the equation of motion 

pu.Au = C {t)^ AR + (j l+ - - l(h(t))) (36) 
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An efficient method to find the noise spectrum is provided by the MacDonald formula, 22 
which has been widely used in quantum transport. 21 ' 23-25 The validity of this formula requires 
that the current correlation function (7(t 1 )/(t 2 )) is only a function of the time difference t\ — 
t 2 and that, therefore, the transport is in the stationary regime. A direct application of the 
MacDonald formula to the present time-dependent transport problem is thus not possible. 
However, in the present study the driving field is time-periodic. The discrete temporal 
translation symmetry H{t + T) = H(t) makes it possible to generalize the MacDonald 
formula for the noise spectra time-averaged over one period as 

SiAu) = 2c 2 f T dt Tldot \ S {- lU j)^ +& ^^ +s ^(-iuj) - S{iu)^ +S ^^ +S ^\iu) \ , (37) 
<jj / Jo 2i L J 

where the composents of the superoperator S(s) are given by 

[S(s)]lmk;l'm'k' = ( s ~ Su'5 mm '5 k k' ■ (38) 

To clarify the meaning of this definition, we note that the superoperator S(s) acts on an 
arbitrary operator A with Fourier-transformed matrix elements Ai m -k in our standard basis 
|0,|m) = |0),|t>,U),m> as 

[iS(s) A]i m -k = ^2 *^2[<S(s)}l m k;l> m >k' M'm'^k'- (39) 
I'm 1 k' 

The derivation of the generalized MacDonald formula is outlined in the appendix. In evalu- 
ating the current noise from the generalized MacDonald formula, we encounter the Laplace 
transforms J rm,n (s) of the moments of the electron number operators in the leads. These 
moments are nothing but the expansion coefficients of J 7 defined in Eq. (15). 

Taking the trace over the dot degrees of freedom and the average over one period makes 
only the matrix elements J-^^q of J :m ^ n that are diagonal in the dot basis and have Floquet 
index k = contribute to the final result. We arrive at the expression 

= _ e 2 w Tr dot [J- fc tV 5i ' L,<5m+<5 " fl (-^) + Fto 5 " 1 -' 6111 ^" 11 ^)}. (40) 

CO 

Now we require the charge moments in the left and right leads. They can be obtained from 
the equation of motion for J 7 [Eq. (19)]. Suppose we switch on the counting fields at some 
time ti, before that time the system can be described by the density matrix without the 
counting fields in the (quasi-) stationary limit. Since we have assumed to - ► — oo above, 
any initial correlation have died out at time ti. 49 We set the number of electrons having 
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tunneled into lead / up to time t± to zero, Ni(ti,ti)=0. After time t±, the system evolves 
under the influence of the counting fields. Then, we solve the equations for J-" m ' n (t) by means 
of Laplace transformation. For example, the solution of Eq. (23) reads 



JF°'°(s) = [S(s)-C]- l ^°(h 



(41) 



where J- ' (ii) can be found from the stationary master equation CJ 70,0 = in the absence 
of counting fields. Analogously, we find expressions for the other expansion coefficients after 
the Laplace transformation as 



jF x >° = {S-C)- 1 {Ji+- Ji-)F 



0.0 



(42) 



0,0 



(43) 



0.0 



(44) 



P>> 2 = (5 - Z:)" 1 [2(^2+ " J2-)J' ' 1 + (J2+ + J2-)F 



0.0 



(45) 



f 1 ' 1 = {s- cy l Uj 1+ - j^)^ 1 + {j 2+ - j 2 -)t 



1.0 



(46) 



where we have omitted the arguments s. The solutions for these coefficients together with 
the generalized MacDonald formula [Eq. (37)] give the desired time-averaged current-noise 
spectra of the ac-driven quantum dot. The approach presented in this study can easily be 
generalized to take more complex structures with multiple levels and inter-level transitions 
into account. 



III. RESULTS AND DISCUSSION 

In the following, we present our numerical results based on the Floquet master-equation 
method and discuss the transport properties of the single-level quantum dot with time- 
dependent fields. An additional dc magnetic field along the z or x direction is taken into 
account; it splits the energy levels of the singly charged quantum dot due to the Zeeman 
effect. In the present study, we choose the ac field to be either a rotating magnetic field 
in the xy plane or an ac gate voltage. The ac gate voltage and the rotating magnetic 
field will affect the quantum conductor in quite different manners. An ac gate voltage 
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only changes the eigenvalues of H^ ot (t) periodically, which in the adiabatic limit of large 
T become the eigenenergies. The ac gate voltage will not induce any transition between 
different eigenstates of H dot (t) (no spin flip is possible) because the eigenstates are unaffected 
by the gate voltage. The electron is trapped in one spin state. The situation is different 
for a rotating magnetic field. A rotating magnetic field does not change the eigenvalues of 
Hdot(t) but does change the eigenstates and thus can flip the spin of the electron. The spin 
polarization of the dot will thus evolve with the rotating magnetic field. The two types of ac 
fields show drastically different behaviors in the transport properties as we will show below. 

The full Hamiltonian of the quantum dot is written as (we reiterate that we choose 
\e\ = h = k B = 1), 

#dot = J2( eV c + Vec cos ttt + aB z ) d\d a + U d\d^d\d± 

+ B x (4dj. + d\d t ) + B ac {d\d^ m + d\d t e~ int ) , (47) 

where cVq is the on-site energy of the quantum dot due to the dc component of the gate 
voltage Vq, Kc is the amplitude of the oscillating gate voltage, U represents the intra-dot 
Coulomb interaction, and B x and B z are half the Zeeman energies of the singly occupied 
dot due to the dc magnetic fields in the x and z direction, respectively. Half the Zeeman 
energy of the rotating magnetic field is given by B ac . Note that while it is customary to talk 
about photon-assisted processes in this context, the treatment of the electromagnetic field 
in the Hamiltonian is completely classical. 

We work in the sequential-tunneling regime and choose a symmetric coupling geometry 
with a = b. We assume that the bias voltage Vd c symmetrically shifts the chemical potentials 
by Hl,r = ±eVdc/2. In the framework of wide-band approximation, the tunneling rate is 
given by r^ n;ij (e) = 27rV^ j(T V^™*, where we have assumed the coupling strength V^ k cr = V 
to be a constant and have set the density of states of lead I to unity. In the present study, 
we have assumed the tunneling matrix to be independent of the energy and the occupation 
number on the dot. An inclusion of state-dependent tunneling is straightforward. We 
assume that the electrons tunneling in and out of the dot with an energy-independent rates 
T = r^ n .^(e) and set T — 1 as the energy unit. 
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A. Differential conductance 



We start our discussion with the differential conductance. The gray-scale plot Fig. 1 
shows the differential conductance dI/dV^ c vs. the dc bias voltage V& c and the gate voltage 
Vq with or without an ac field. The calculations are for the Coulomb interaction strength 
U = 24 and the temperature ksT = 0.32. The frequency of the ac field is Q = 8. 

Without an ac field, Fig. 1(a) gives the familiar diamond structure due to the Coulomb 
blockade. Numerical results for the differential conductance when the quantum dot is mod- 
ulated by an ac gate voltage are presented in Fig. 1(b). Fig. 1(c) gives the results when 
the quantum dot is modulated by a rotating magnetic field. Fig. 1(d) shows the differential 
conductance when the quantum dot is modulated by a rotating magnetic field in the xy 
plane while a dc magnetic field is applied in the x direction, i.e., in the plane of the rotating 
magnetic field. 

When there is an ac field, several striking features emerge in the differential conductance: 

(1) At the edge of the Coulomb diamond, the sharp differential-conductance peak for the dc 
transport shown in Fig. 1(a) is partially suppressed by the ac gate voltage or the rotating 
magnetic field. Note the different gray scales in Fig. 1 (a), (b), (c), and (d). This can be 
attributed to the suppression of the elastic resonant peak by the photon-assisted processes. 

(2) In the presence of an ac field, there are lines parallel to the edges of the Coulomb diamond. 
The distance of these lines to the peak position is approximately the frequency of the ac 
field, indicating a photon-assisted tunneling process. (3) An interesting feature of these 
lines can be observed inside the Coulomb diamonds. For the ac gate voltage, the Floquet 
quasienergies are spin degenerate. Therefore, the main lines in the differential conductance 
plot in Fig. 1(b) are not split. However, satellites due to photon-assisted inelastic tunneling 
events appear, in which an energy quantum of Q is absorbed from or emitted into the driving 
field. We see from Fig. 1(b) that these additional lines remain distinct inside the Coulomb 
diamond. On the other hand, when the quantum dot is driven by a rotating magnetic field, 
the quasienergies are not degenerate. Therefore, the main elastic lines are split into two at 
the edge of the Coulomb diamond in Fig. 1(c). Interestingly, the lines due to the photon- 
assisted tunneling now only appear outside of the Coulomb diamond, as can be seen in Fig. 
1(c), indicating that the photon-assisted tunneling is forbidden inside the Coulomb diamond. 
When the quantum dot is modulated by a rotating magnetic field and a dc magnetic field is 
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applied in the plane of the rotating field, these lines in the Coulomb-blockade regime revive. 
This can be clearly seen in Fig. 1(d). 

The disappearance of the photon-assisted tunneling inside the Coulomb diamond for a 
pure rotating magnetic field can be understood as follows. In the Coulomb diamond, the 
Floquet quasienergies corresponding to the singly occupied states are far below the Fermi 
energies of the two leads. A direct tunneling between the dot and the leads is forbidden due 
to the Pauli principle and Coulomb blockade. Therefore, an electron is effectively trapped in 
one quantum state on the dot. According to our previous discussion, only the spin direction 
of this quantum state can evolve with the rotating magnetic field. However, its eigenvalues 
of Hdot{t) remain unchanged. Therefore, the electron cannot gain extra energy from the ac 
magnetic field. As a consequence, we cannot observe lines due to photon-assisted tunneling 
inside the Coulomb diamond. Outside of the Coulomb diamond, the tunneling between the 
dot and the leads becomes possible. When an electron is injected from the lead into the dot, 
the system can absorb or emit photons, i.e., the Floquet index k can change. One could say 
that transport happens via several Floquet channels. Such photon-mediated tunneling can 
then give rise to the photon-assisted differential-conductance peaks. 

The situation becomes different when a dc magnetic field is applied in the plane of 
the rotating magnetic field as shown in Fig. 1(d). In that case, the eigenstates and the 
eigenvalues of H^ ot (t) change periodically. Electrons can gain extra energy from the ac 
field by absorbing or emitting a photon. In the Coulomb diamond, electrons on the dot 
are able to tunnel out via the photon-assisted tunneling and we again find the lines due to 
the photon-assisted differential-conductance peaks inside the Coulomb diamond as shown in 
Fig. 1(d). 

B. Zero-frequency Fano factor 

In the following, we show numerical results for the time-averaged zero-frequency noise 
of the quantum dot. The zero- frequency noise has been studied by the quantum master- 
equation method in the stationary 12 ' 16,23 and also in the time-dependent case. 44 Without 
ac field and at zero temperature, the zero-frequency Fano factor S(0)/2eI describes the 
deviation of the shot noise from its Poissonian value. We choose the parameters T = 0.32, 
U = 8, eVc = 8, and Q = 4.8. We assume that a dc magnetic field in the z direction, 
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B z = 1.6, is applied to the quantum dot. The finite value of U and the Zeeman splitting 
make it possible to see plateaus in the Fano factor at different occupation numbers on the 
dot. 16 Fig. 2 shows the zero-frequency Fano factor as a function of the dc bias with or without 
an ac gate voltage. 

Without an ac gate voltage, V ac = 0, the results reproduce the main features reported 
in Ref. 16. At very low bias voltage Va c — > 0, the main contribution to the noise is the 
finite thermal noise while the current as well as the shot noise are suppressed. Therefore, 
the Fano factor diverges at V^ c = 0. For low dc bias voltage, the quantum dot operates in 
the Coulomb-blockade regime. With further increasing dc bias voltage, the energy levels of 
the quantum dot one by one enter the transport window defined by the dc bias. This can 
be clearly identified in the Fano factor by the plateaus at different values. The edges of 
the plateaus are broadened by the finite temperature. For very large dc bias, where all the 
energy levels of the quantum dot lie in the transport window, the Fano factor approaches 
the well-known limit of 1/2 for our symmetric-coupling case. 16 

The results for the dc case demonstrate that the plateaus of the Fano factor can give a 
good account of the transport channels. 16 In the presence of an ac field, we now consider 
the time-averaged Fano factor 5(0)/2e7. We can see from Fig. 2 that for small Vd c the Fano 
factor becomes larger as we increase the amplitude of the ac gate voltage. On the other hand, 
additional photon-assisted transport channels are available due to the ac field, which will 
modify the Fano-factor curve. With increasing ac field, the Fano factor will thus deviate from 
the plateau behavior seen in dc case due to the opening of these photon-assisted transport 
channels. When the ac gate voltage is large enough, the Floquet eigenstates that lie outside 
of the transport window can contribute to the current via photon-assisted tunneling. As 
a consequence, the plateaus in the Fano-factor curve become vague. For very large bias 
voltages, all the Floquet levels are well inside the transport window. The Fano factor then 
will approach the same value 1/2 as for the time-independent transport. In Fig. 3, we present 
our results for the zero-frequency current noise in the presence of a rotating magnetic field. 
As in the case of an ac gate voltage, the Fano factor deviates from the dc behavior with 
increasing ac field. Additional plateaus can be observed in the Fano-factor curve when 
we vary the dc bias voltage. Transitions between plateaus result from additional Floquet 
channels becoming available. 
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C. Frequency-dependent Fano factor 

Now we present our results for the full current-noise spectra of a quantum dot under an 
ac field. The noise spectra have previously been studied in the stationary-transport regime. 
An analytical expression for the noise spectrum of a single-level quantum dot can be found 
in Ref. 23. Unless stated otherwise, the following calculations assume U = 24, V& c = 12.7, 
T = 1.6, Q = 8 and eVc = —8. We introduce the frequency- dependent Fano factor S(u)/2el 
to characterize the time-averaged noise power. As discussed previously, the ac gate voltage 
and the rotating magnetic field will modulate the quantum conductor in different ways. In 
the following, we show that the noise spectra are also strikingly different. 

In Fig. 4, we present the results for the frequency-dependent Fano factor S(u)/2el as a 
function of the frequency uj for different amplitudes V ac . No dc magnetic field is applied. 
Without an ac field (V ac = 0), the noise spectrum shows a peak at zero frequency and 
approaches a constant value for large u. The peak in the noise spectrum is due to the elastic 
processes in the transport. 29 When an ac gate voltage is applied, additional structures in 
the noise spectra are expected due to photon-assisted processes. For the present set of 
parameters, one can clearly see that with increasing amplitude of the ac gate voltage, an 
additional peak appears in the noise spectrum. While the height and width of this peak 
vary a lot with increasing amplitude, its peak position u p remains almost unchanged at the 
external driving frequency Q. 

Now we turn to the rotating magnetic field in the xy plane. In Fig. 5, we plot the Fano 
factor as a function of the frequency u for different amplitudes B ac of the rotating magnetic 
field. Similarly to the results presented in Fig. 4, a peak is generated and the width and 
height of this peak depend on the amplitude. However, the peak position is not fixed at Q 
in contrast to what we have observed in Fig. 4 for the ac gate voltage. Instead, its position 
shifts with increasing amplitude, as shown in Fig. 5. 

By comparing Fig. 4 and Fig. 5, we see that the peak position of the noise spectra behaves 
differently when we increase the ac strength, depending on the type of the ac field. Recalling 
that when electrons tunnel through a time-independent quantum two level system, its current 
noise spectra show additional structure at the energy difference of the two transport channels 
of the system due to its internal coherent dynamics, 14 we will show that the peak position 
of the noise spectra for ac transport can be understood from the interference between two 
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possible Floquet transport channels. If the quantum dot is modulated by a rotating magnetic 
field, the last term in the dot Hamiltonian [Eq. (47)] shows that the ac magnetic field 
couples one spin state with the quasienergy e with a state with the opposite spin and the 
quasienergy e — Q (in the extended zone scheme). The coupling strength is given by _B ac . 
The corresponding Floquet Hamiltonian then decomposes into 2x2 blocks of the form 

h Fl =( 6 BaC ) . (48) 
\B ac e-n) 

The resulting quasienergies in the first Brillouin zone [0, Q) are 



ei = e~2+ 1 2 ' ( } 



3 fi ^ + ABl 



with the difference 



co P = e 2 - ex = 20. - JQ? + 45 a 2 c (51) 



(these expressions hold if B ac <y/3Q/2). 

If now an electron tunnels into the dot, the system ends up in a superposition of the two 
Floquet states, the phases of which change with different angular frequencies, corresponding 
to spin precession with the difference frequency u p . When the electron tunnels out again, 
the superposition is projected onto the spin direction of the original electron since lead 
electron creation and annihilation operators are paired with identical quantum numbers in 
the master equation. This leads to interference with a typical frequency u p , which enhances 
the current- current correlation function Sij(t, if) in Eq. (30) for t — t' being a multiple of the 
period 2tt/uj p and thus leads to a peak in the noise spectrum at u p . The peaks seen in Fig. 
5 are indeed centered at ui p given by Eq. (51). 

Comparing with the stationary transport through a stationary two level system, 14 the 
transport through a quantum dot with rotating magnetic field can be understood as another 
type of two level quantum system. The significant difference here is that our two levels are 
defined by the Floquet channels due to a periodic ac field and not by the true eigenenergies 
of an time-independent Hamiltonian. 

When an ac gate voltage is applied to the quantum dot, the ac field will not couple the 
different spin states. Only the eigenvalues of H dot (t) will be modulated, see Eq. (47). The 
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corresponding Floquet Hamiltonian decomposes into two infinite blocks for the two spin 
directions, where each block has the tridiagonal form 



\ 



e + tt Kc/2 
Kc/2 e Kc/2 
Kc/2 e - n 



\ 



(52) 



Therefore, the electrons can tunnel through the quantum dot via infinitely many Floquet 
channels with the same quasienergy in the first Brillouin zone [0, Q) but all possible Floquet 
indices k. The quasienergies in the extended zone scheme thus differ by integer multiples 
of Q. These quasienergy differences define the peak positions in the noise spectra. In Fig. 
4, a peak at Q appears, corresponding to two channels with their Floquet indices (photon 
numbers) differing by unity. One should also expect peak structures at nQ, n > 1. However, 
to observe these peak structures, one may need a stronger ac field to enable multi-photon- 
assisted transport. In the inset of Fig. 4, a small shoulder emerges at 2Q for the largest 
amplitude, Kc = 8. 

So far in our discussion, no dc magnetic field has been considered. For the quantum dot 
with an ac gate voltage and a dc magnetic field in the z direction, our results of the noise 
spectra for different voltage amplitudes are displayed in Fig. 6. The parameters are the same 
as those used in Fig. 4 except that the strength of the dc magnetic field in the z direction 
is B z = 1.6. For the present set of parameters, the peak at Q is replaced by a dip. We 
observe that the appearance of the peak or the dip depends on the detailed parameters used 
in our calculation. The peak or dip position remains unchanged as we increase the voltage 
amplitude. We have checked that the noise spectrum does not depend on the direction of 
the dc magnetic field. This is because the full SU(2) symmetry is preserved since we have 
included the time evolution of the off-diagonal elements of the reduced density matrix within 
our quantum master-equation approach. 

Contrary to the quantum dot with an ac gate voltage, for the case of a rotating magnetic 
field, the noise spectra do depend on the direction of the additional dc magnetic field. When 
the dc magnetic field is perpendicular to the plane of the rotating magnetic field, the noise 
spectra behave much like those for a pure rotating magnetic field. Only one peak appears 
at a non-zero frequency and the peak position shifts with the amplitude of the rotating 
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magnetic field. Numerical results for the noise with a rotating magnetic in the xy plane and 
a dc magnetic field in the z direction are displayed in Fig. 7. The parameters are the same 
as in Fig. 5 and the dc magnetic field is B z = 1.6. It is easy to verify that the peak position 
can again be determined by the difference between two Floquet quasienergies. 

When a non-zero dc magnetic field B x is applied in the plane of the rotating magnetic 
field, the Floquet Hamiltonian cannot be reduced to a 2 x 2 matrix form as for the previously 
discussed situation of vanishing dc magnetic field, since the B x term in Eq. (47) mixes spin- 
up and spin-down states. Together with the rotating magnetic field this couples all Floquet 
states with the same quasienergy in the first Brillouin zone and different Floquet indices. 
Thus the electrons can tunnel through the quantum dot via infinitely many Floquet channels. 
The interference between these Floquet channels then gives rise to much richer behavior in 
the noise spectra. Numerical results for the noise spectra of a quantum dot driven by a 
rotating magnetic field in the xy plane and with a dc magnetic field in the x direction are 
presented in Fig. 8. The parameters used in the calculation are the same as in Fig. 5 and 
the dc magnetic field in B x = 1.6. In Fig. 8, more peaks are observed than for vanishing 
dc magnetic field in Fig. 5. One can see that besides the peak determined by Eq. (51), 
there are both peaks (dips) fixed at Q and structures at positions depending on the ac-field 
amplitude B ac . All peak (dip) positions correspond to the differences of available Floquet 
quasienergies. 



IV. SUMMARY 



In this paper, the transport properties of a single-level quantum dot modulated by ei- 
ther an ac gate voltage or a rotating magnetic field have been studied within the Floquet 
quantum master-equation approach in the sequential-tunneling limit. We have employed a 
generalized MacDonald formula to obtain the time-averaged current noise spectra for both 
cases. Numerical results for the differential conductance and the frequency-dependent cur- 
rent noise have been presented. Besides the usual diamond structure due to the Coulomb 
blockade in the differential conductance, photon-assisted tunneling can give rise to addi- 
tional lines parallel to the edges of the Coulomb diamond. These lines cannot survive inside 
the Coulomb diamond in the case of a rotating magnetic field. This is due to the fact that 
the rotating magnetic field only periodically rotates the spin direction while the energy of 
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the electron on the dot remains unchanged. The frequency- dependent noise spectra of the 
quantum dot show additional peaks or dips in the presence of an ac field. The behavior 
of these additional structures depends on the nature of the ac driving field. In the case 
of an ac gate voltage, the position of the finite-frequency peak is fixed at the external ac 
frequency, independently of the voltage amplitude. On the other hand, in the case of a 
rotating magnetic field, the peak at non-zero frequency moves with changing amplitude of 
the rotating magnetic field. An additional dc magnetic field in the plane of the rotating 
magnetic field can also drastically change the noise spectra: it leads to the appearance of 
both movable and fixed peak structures in the noise spectra. All these peak positions are 
found to be determined by the energy differences between two Floquet transport channels. 
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Appendix A: Derivation of the generalized MacDonald formula with ac field 



A derivation of the MacDonald formula for time-independent transport has been pre- 
sented in Ref. 24. We show here that the periodicity of our time- dependent Hamiltonian 
makes it possible to estimate the noise spectra by generalizing the MacDonald formula. 
The derivation of the generalized MacDonald formula for the time-averaged noise spectra is 
outlined in the following. 

We start from the Fourier-transformed current correlation function 



S{t,uj) 



dre WT ({5I(t),5I(t-T)}} 



(Al) 



where 5I(t) = I(t) — (I{t)) and {A, B} = AB + BA. For simplicity, we omit the lead index 
in the noise expressions in this appendix. 
From the definition of the current, we have 
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where N(t + r, t) denotes the number of charges transferred during the interval from t to 
t + r. Taking the expectation value of the square of this equation, we obtain 
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Inserting the Fourier decomposition of the time- dependent relation 
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into Eq. (A3), we obtain 
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where we have used the notation J2'k = J2k^o- Differentiation with respect to r gives 
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We next perform a Fourier transformation and take the time average over one period. 
The second term on the right-hand side vanishes due to its periodicity. Since the cur- 
rent correlation function is symmetric, S(t,t') = S(t',t), it can be shown that S(u) = 
T lo dt j'^ 00 dt'S{t,t')e iu ' { - t ~ tl) has the property S(ui) = S(-u). We arrive at the generalized 
formula for the time-averaged noise spectrum for a periodic driving field, 
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Noting that the integrand of the r integral is even, we obtain the final result for the gener- 
alized MacDonald formula for the time-averaged noise spectrum, 
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In comparison to the MacDonald formula for steady state transport, an integration of t over 
one period is carried out to obtain the time-averaged noise spectra. The time average can 
also be expressed by an average over the initial phase of the ac field. Hence, our expression 
is equivalent to the form given by Clerk and Girvin. 46 
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FIG. 1. Gray-scale plots of the differential conductance dl/dV^c as a function of the dc bias voltage 
and the dc gate voltage for a quantum dot (a) without any ac fields, (b) modulated by an ac gate 
voltage, (c) modulated by a rotating magnetic field, and (d) modulated by a rotating magnetic field 
and with an additional dc magnetic field applied in the plane of the rotating field. Dark regions 
represent low differential conductance. The frequency of the ac field is f2 = 8. The amplitude of 
the ac gate voltage in (b) and of the rotating magnetic field in (c), (d) are V ac = 6.4 and B ac = 3.2, 
respectively. The dc magnetic field in (d) is B x = 1.6. 
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FIG. 2. (Color online) Fano factor of the quantum dot as a function of dc bias voltage Vd c in unit 
of Vq for different amplitudes V &c of the ac gate voltage. The parameters of the device are given 
in the text. 
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FIG. 3. (Color online) Fano factor of the quantum dot as a function of dc bias voltage V^ c in unit 
of Vg for different amplitudes B &c of the rotating magnetic field. The other parameters are the 
same with those of Fig. 2. 
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FIG. 4. (Color online) Noise spectra for different gate voltage amplitudes V ac . Independently of 
V ac , the main peak position is fixed at fi. The inset shows an enlarged view of the noise spectra 
around 2Q . 



32 




frequency co/Q 



FIG. 5. (Color online) Noise spectra with a rotating magnetic field in the xy plane with various 
amplitudes B &c . The peak position shifts with increasing B ac . 
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FIG. 6. (Color online) Noise spectra for a quantum dot modulated by an oscillating gate voltage 
with various amplitudes V ac in the presence of a dc magnetic field B z = 1.6 in the z direction. The 
dip position is fixed at the driving frequency. 
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FIG. 7. (Color online) Noise spectra in the presence of a rotating magnetic field in the xy plane 
with various amplitudes B ac . A dc magnetic field B z = 1.6 is applied in the z direction. Only one 
peak appears, which shifts with B ac . 
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FIG. 8. (Color online) Noise spectra with a rotating magnetic field in the xy plane with various 
amplitudes B ac . A dc magnetic field B x = 1.6 is applied in the x direction. New features appear 
in the noise spectra since more Floquet channels are involved in the transport in the presence of 
an in-plane magnetic field. 
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